o 

(N 



O 



The Astrophysics of Planetary Systems: Formation, Structure, and 
Dynamical Evolution 

Proceedings IAU Symposium No. 276, 2010 © 2010 International Astronomical Union 

A.C. Editor, B.D. Editor & C.E. Editor, eds. DOI: 00.0000/X000000000000000X 

Non-convergence of the critical cooling 
timescale for fragmentation of 
self-gravitating discs 

Farzana Meru 1,2 f and Matthew R. Bate 1 

1 School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL, UK 

2 Institut fur Astronomie und Astrophysik, Universitat Tubingen, Auf der Morgenstelle 10, 

72076 Tubingen, Germany 



Abstract. We carry out a resolution study on the fragmentation boundary of self-gravitating 
discs. We perform three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations of 
discs to determine whether the critical value of the cooling timescale in units of the orbital 
timescale, /3 c rit, converges with increasing resolution. Using particle numbers ranging from 31,250 
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1* I , to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, 
' fragmentation occurs for longer cooling timescales as the resolution is increased. These results 
I certainly suggest that /3 C rit is larger than previously thought. However, the absence of conver- 
gence also questions whether or not a critical value exists. In light of these results, we caution 
against using cooling timescale or gravitational stress arguments to deduce whether gravitational 
instability may or may not have been the formation mechanism for observed planetary systems. 
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■ 1. Introduction 

There are two quantities that have historically been used to determine whether a self- 
gravitating disc is likely to fragment. The first requires the stability parameter, Q < 1 
t— I . (Toomre 1964), where Q — , c s is the sound speed in the d the epicyclic 

frequency, which for Keplerian discs is rs fi, the angular frequency, £ is the surface mass 
density and G is the gravitational constant. Gammie (2001) showed the need for fast 
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\T) ■ u(^^) O, where u is the specific internal energy and du coo i/di is the total cooling 



cooling and suggested that if the cooling timescale can be parameterised as (3 — t coo \Q 
u {~d± 2L ^) wnere u is the specific internal energy and du coo i/di is the total coolir 
rate, then fragmentation requires (3 < /3 C rit, where the critical cooling timescale, /5 cr it ~ 3 



for a ratio of specific heats, 7 = 2. Rice et al. (2005) performed 3D SPH simulations to 



show the dependence of /3 C rit on 7: for discs with 7 = 5/3 and 7/5, /3 C rit ~ 6 — 7 and 
w 12 — 13, respectively. The above authors show that the cooling condition is equivalent 
to a maximum gravitational stress that a disc can support without fragmenting. Recently, 



Meru & Bate (2010) suggested that /3 cr i t depends on the disc and star conditions 



We carry out a thorough convergence test of the value of /3 cr it required for fragmenta- 



tion using a 3D SPH code. We simulate a 0.1M© disc surrounding a IMq star, spanning 
a radial range, 0.25 < R ^ 25 au, using 7 = 5/3. The initial surface mass density and 
temperature profiles are £ tx Rr 1 and T tx i? -1 / 2 , respectively, and initially, Q > 2 ev- 
erywhere. We carry out simulations with this disc setup using 31,250, 250,000, 2 million 
and 16 million particles. We simulate the discs using various values of f3 to determine 
/3 cr it at different resolutions. 
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Figure 1. a. Surface mass density rendered images of the fragmenting discs with 31,250, 250,000, 
2 million and 16 million particles. At higher resolution, the disc can fragment for larger /3. The 
axes scale from -25 au to 25 au in both directions, b. Graph of /3 against resolution of the 
non- fragmenting (open squares), fragmenting (solid triangles) and borderline (open circles) sim- 
ulations (defined as discs that fragment but quickly shear apart with no further fragmentation). 
The solid black line divides the fragmenting and non-fragmenting cases and the grey region is 
where fragmentation can occur. The graph shows no evidence of convergence with resolution. 
The thin dotted line shows how the trend will continue if convergence is not reached with even 
higher resolution. If convergence can be achieved, the dotted line would follow a flatter profile. 

2. Results 

Figure shows images of the fragmenting discs at various resolutions. We see that 
fragmentation occurs for higher values of j3 as the resolution increases. Figure Q}} sum- 
marises all the simulations performed. With the data that is available, the dividing line 
between the fragmenting and non fragmenting cases increases linearly with linear reso- 
lution and therefore convergence has not been reached. 



3. Numerical and observational implications 

The lack of convergence certainly shows that the critical value of the cooling timescale 
is, at the very least, longer than previously thought. However, it also opens up the 
possibility that there may be no value of /3 for which such a disc can avoid fragmentation, 
given sufficient resolution. If this is the case, it suggests that the problem may be ill posed. 
In other words, it may not be possible for a disc to settle into an equilibrium where there 
is a balance between heating from gravitational instabilities and a simple imposed cooling 
timescale. This implies that a self-gravitating disc that cools at a rate given by = 
and is only heated by internal dissipation due to gravitational instabilities, may not be 
able to attain a self-regulated state and will always fragment, regardless of the value of 
p. This re-opens the question of what the criterion for fragmentation of a self-gravitating 
disc really is and in addition, where in a disc fragmentation can realistically occur. These 
results cast some serious doubts on previous conclusions concerning fragmentation of self- 
gravitating discs. In addition, since cooling timescale arguments can be used to determine 
at what radii in a disc fragmentation can occur (e.g. Rafikov 2009; Clarke 2009), the 
results presented here need to be considered when making conclusions as to whether 
observed planetary systems may or may not have formed by gravitational instability. 




Resolution effects on fragmentation 
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